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We carry out a theoretical investigation to clarify the dynamic property of photo-created 
nano-sized ferroelectric cluster observed in the paraelectric BaTiOa as a real time correla- 
tion of speckle pattern between two soft X-ray laser pulses, at just above the paraelectric- 
ferroelectric phase transition temperature. Based on a model with coupled soft X-ray 
photon and ferroelectric phonon mode, we study the time dependence of scattering 
probability by using a perturbative expansion approach. The cluster-associated phonon 
softening as well as central peak effects are well reproduced in the phonon spectral func- 
tion via quantum Monte Carlo simulation. Besides, it is found that the time dependence 
of speckle correlation is determined by the relaxation dynamics of ferroelectric clus- 
ters. Near the transition point, cluster excitation is stable, leading to a long relaxation 
time. While, at high temperature, cluster structure is subject to the thermal fluctuation, 
ending up with a short relaxation time. 

Keywords: Nonlinear optics; soft X-ray laser; nano cluster. 



1. Introduction 

Ferroelectricity refers to the ability of developing a spontaneous polarization (elec- 
tric dipole moment) on decreasing temperature below Curie point (T c ). As a proto- 
type of the ferroelectric perovskite compounds, barium titanate (BaTiOa) under- 
goes a transition from paraelectric cubic to ferroelectric tetragonal phase at T c =395 
K. In addition to its extensive application in the technology such as non- volatile 
memory devices, capacitor applications and piezoelectric transducers due to its high 
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dielectric constant and switchable spontaneous polarization^ there is also an en- 
during interest in understanding the mechanism of paraelectric-ferroelectric phase 
transition in it. Although it is generally considered that the transition might be 
a classic displacive soft-mode type driven by the anharmonic lattice dynamics 1213J 
recent studies have suggeste d an existence of order-disorder instability competing 
with the displacive transitiorP^, making this issue still controversial to date. 

Direct observation on creation and evolution of ferroelectric cluster around T c 
is of crucial importance for clarifying the nature of phase transition. However, 
conventional techniques, like neutron scattering,^ polarizing optical microscopy-^ 
and scanning probe microscopyp 1 can hardly detect the ultrafast transient status of 
dipole clusters. In contrast, diffraction speckle patterns of BaTiC>3 crystal measured 
by the soft X-ray laser has turned out to be an efficient way to elucidate the 
dynamics of domain structure. 9-11 

Very recently, NamikawJ^ study the nano-sized polarization clusters in BaTi03 
at just above T c by the plasma-based X-ray laser speckle measurement. In this 
experiment, two consecutive coherent X-ray laser pulses of 160 A with an adjustable 
time difference are generated by the Michelson type beam splitter. After the photo 
excitation by the first pulse, ferroelectric clusters of nano scale are created in the 
paraelectric BaTiC>3 and tends to be smeared out gradually due to the thermal 
fluctuation. This relaxation of cluster structure thus can be recorded in the speckle 
patterns of the second pulse as a function of its delay time from the first pulse, 
which is known as the speckle intensity correlation function. It has been found that 
near the T c , this correlation follows an exponential decay process, i.e., the intensity 
of speckle pattern declines as the delay time increases. Moreover, the decay rate 
also decreases upon approaching T c , indicating a critical slowing down of the dipole 
relaxation time. Hence, by measuring the correlation between two soft X-ray laser 
pulses, the real time relaxation dynamics of polarization clusters in BaTiC>3 is 
clearly represented. 

In this work, we examine the above-mentioned novel behaviors of ferroelectric 
cluster observed by Namikawa from a theoretical point of view, aiming to provide 
a basis for understanding the critical nature of BaTiC>3. Our calculation confirms 
that the relaxation of photo-created nano clusters plays an essential role in deter- 
mining the delay time dependence of speckle correlation function. The remaining 
of present paper is organized as follows. In Sec. 2, the model Hamiltonian and the- 
oretical treatment are put forward. In Sec. 3, our numerical results are discussed in 
comparison with the experimental discoveries. In Sec. 4, a summary with conclusion 
is presented finally. 



2. Theoretical model and method 

In order to describe the optical response of BaTiC>3 due to X-ray scattering, we 
design a theoretical model to incorporate the radiation field and ferroelectric dipole 
correlation as well as a weak interplay between them. As mentioned above, the phase 
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transition in BaTiOa is likely to be an admixture of displacive and order-di sord er 
types. For this sake, we adopt a Krumhansl-Schrieffer type phonon mode^l to 
characterize the ferroelectric anharmonicity in BaTiOs, for this model is suitable for 
describing the crossover between displacive and order-disorder phase transitions!^ 
In this scenario, our model Hamiltonian reads (h — 1 and ke = 1 throughout this 
paper): 

H = Hp + Hf + H p f, 

Hp = ^k%ak, 
k 

i ^ ^ l ' <i,v> 

V \ - + 

H Pf = N Z-^ a fc + |«fe-f Qq'-iQ-q'-i , (2- 1 ) 
q,q',k 

where at (a*;) is the creation (annihilation) operator of a photon with a wave 
vector k and an energy Qk- Qi is the coordinate operator for ionic displacement 
along the easy axis of BaTiC>3 at the site I with a dipole oscillatory frequency 
ujq. Q q (= N^ 1 / 2 e~ lql Qi) is its Fourier component for wave vector q. Without 
losing generalitivity, here we use a simple cubic lattice, and the total number of 
lattice site is N. Hf is the extended Krumhansl-Schrieffer model for phonon, where 
a sixth order term of Qi is introduced to produce a triple- well potential, because the 
phase transition of BaTiOa is first order. Thus, the lattice dynamic properties are 
determined by the frequency u>q, temperature T, and dimcnsionless parameters C4, 
cq, d%. Hpf corresponds to a Raman type bi-linear coupling between the ferroelectric 
phonon mode and photons as the phonon mode has an even parity in the present 
case, and V is the coupling strength. 

Since there are two photons involved in the scattering with crystal, the photon- 
phonon scattering probability can be written as, 

P(t) - 2((a fc0 (0)oi 1 (A)a fc0 (t)a^(A + t) 

xa fc ,(A + t)4 o (t)a fel (A)4 o (0))), (2.2) 

where, 

((•••)) =Tr(e-^-..)/(e-^), (2.3) 

means the expectation, (3 (= 1/T) is the inverse temperature, and the time depen- 
dent operator 0(t) is defined in the Heisenberg representation, 

0(t) = e ltH Oe- ltH . (2.4) 

Here, t denotes the time difference between two incident X-ray laser photons, and 
fco is the wave vector of them. After a small time interval A, these photons are 
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scattered out, respectively, k\ and k[ are the wave vectors of the first and second 
outgoing photon. 

Regarding H p f as a perturbation, we separate Hamiltonian of Eq. (2.1) as, 

H = Ha + H P f , (2.5) 

where 

H = H P +H f , (2.6) 

is treated as the unperturbed Hamiltonian. By expanding the time evolution oper- 
ator in Eq. (2.4) with respect to H p f, 

(2.7) 



e -UH _^ e -itH 



I- if drH pf (r) 
Jo 



we find that the lowest order terms which directly depend on t are of fourth order, 

r A 



P(t)-/ dn f dr 2 f dr[ f dT^((a ko H pf (rl 
Jo Jo Jo Jo ki k , 



xa fcl e- iA ^ J ff f , / (ri)4 o )), ( 2 -8) 

where the operators with carets are defined in the interaction representation, 

O(r) = e iTHo Oe- lTH ". (2.9) 

In Fig. 1, we present a diagrammatic analysis for the phonon-coupled scattering 
process corresponding to Eq. (2.8). In each diagram, the photons (phonons) are 
depicted by the wavy (dashed) lines, and the upper (lower) horizontal time lines 
stand for the bra (ket) vectorsffSI Diagram (a) illustrates the changes of photon 
energy and wave number due to phonon emission and absorption, which corresponds 
to the well-known Stokes and anti-Stokes Raman scattering. Whereas, diagrams (b)- 
(e) represent the exchange, side band, rapid dumping and rapid exchange effects, 
respectively. 

Obviously, diagram (c) bears no i-dependence, while diagrams (d) and (e) only 
contribute a substantial reduction to the time correlation between two laser pulses 
because of the duality of phonon interchange. Accordingly, the i-dependence is 
mainly determined by the diagrams (a) and (b). Allowing for this, the probability 
of interest can be written in the form as, 

pA pA pA pA 

P{t) = ^ d Tl dr 2 dr[ ^_^((a fe0 e^^4 o a fco _ 9 
xe- i(T ^ A) ^at o _ ?a&0 e^^4 o a fco+g e- < ^- A )^4 o+g a feo+9 

t(r 2 -A)H t e -iTjH, a t i(n-A)ff p t e -iTiH p 

AL "fc + 9 " fc o e u fc u fco-9 e a k -q ak o ti 

X 4 )) ((Qg' «)Q q - q > (r{)Q- q + q > (t + T' 2 )Q- ql (t + T' 2 )Q q , (t + T 2 ) 

^ q -At + n)Q-^A-ri)^-An)))- (2.10) 
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Fig. 1. Diagrammatic analysis for scattering process of photon with electric dipole moment 
(phonon). The photons and phonons are denoted by the wavy and dashed lines, respectively. 
In each diagram, the upper and lower horizontal time lines represent the bra and ket vectors, 
respectively. 



Here, we note that the photon and phonon parts have been decoupled, and it can 
been easily seen that the origin of the t-dependence is nothing but the phonon 
(dipole) correlation. 

In the case of forward X-ray scattering, we have fco~fci~fci, thus the normalized 
probability can be further simplified as, 

Pit) _ E q ^\((Q 2 q ))G q+ At)\ 2 

P(o) E q , q >\((Q 2 q Q 2 q+q >W ' 1 ■ ' 

where 



G q (t) = -2i6{t){{Q q {t)Q- q m) - 2»0(-t)«Q_,(O)Q,(t)», (2.12) 
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is the real time correlation function, or real time Green's function of phonon, and 
9(t) is the step function. The Fourier component of this Green's function, 



The phonon spectral function describes the response of lattice to the external per- 
turbation, yielding profound information about dynamic properties of the crystal 
under investigation. Once we get the spectral function, the scattering probability 
and correlation function can be readily evaluated. 

3. Numerical result and discussion 

In the numerical calculation, we set phonon frequency uo=20 meV, and make the 
parameters c 4 =0.076, c 6 =0.00495, d 2 =0.0967. As indicated in Eq. (2.10), the In- 
dependence of probability is a result of phonon correlation. To reveal its relation 
with ferroelectric cluster dynamics, we start our discussion with an inspection on 
the property of phonon spectral function at the temperature regime T>T C . To t his 
end, we first calculate the imaginary time Matsubara Green's function of phonorP21 
on a 12x12x12 cubic lattice with a periodic boundary condition. This is done by a 
path integral Monte Carlo method!^ After obtaining the Matsu bara fu nction, we 
can derive phonon spectral function by the analytic continuation! ^^ 

In Fig. 2, we plot the phonon spectral function in the paraelectric phase at dif- 
ferent temperatures: (a) T=T C , (b) T=1.048T C , (c) T=1.143T C and (d) T=1.429T C . 
In each panel, the spectra are arranged with wave vectors along the TXMR di- 
rection of Brillouin zone (see in the inset), and here u> refers to energy. Since the 
spectra are symmetric with respect to the origin w=0, in this graph we only show 
the positive part of them. When the temperature decreases toward T c , as already 
well-known for the displacive type phase transition, the energy of phonon peak is 
gradually softened. In addition, a so-called central peak, corresponding to the low 
energy excitation of ferroelectric cluster, appears at the T point. On decreasing 
temperature, its intensity is dramatically enhanced, indicating that a stable dipole 
correlation of long-wave-length is built up. As mentioned above, the phonon spec- 
tral function can be understood as the optical response of dipole to the external 
perturbation. Therefore the collective excitation represented by the sharp resonant 
peak in spectra is nothing but the photon-created ferroelectric cluster. 

The appearance of sharp peak near T point signifies that after excitation by soft 
X-ray, the photon-created ferroelectric clusters can be well preserved for a long re- 
laxation time. Thus, they are likely to be observed by subsequent detective photons, 




(2.13) 



is related to the phonon spectral function [= A g (u;)] through jl^l 




(2.14) 
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Fig. 2. Phonon spectral function along the line TXMR in the Brillouin zone at different temper- 
atures of paraelectric phase: (a) T=T C , (b) T=1.048T C , (c) T=1.143T C and (d) T=1.429T C . The 
inset of panel (a) shows the Brillouin zone with high symmetry lines. 
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Fig. 3. Normalized speckle scattering probability as a function of time for paraelectric BaTiOa, 
at T=1.024T C and T=1.048T C . 



resulting in a high intensity of speckle correlation. Keeping this in mind, we move on 
to the results of scattering probability. In Fig. 3, we show the normalized probability 
as a function of time at temperatures T=1.024T C and T=1.048T C . From this figure, 
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we see in both cases, the probability declines exponentially with time, demonstrat- 
ing that the photon-created clusters quickly disappear after the excitation and the 
system relaxes back to the paraelectric state. Meanwhile, at temperature closer to 
T c , more gentle decay rate is observed, indicative of a critical slowing down of the 
relaxation time. This is because with the decrease of temperature, the intensity 
of local polarization is enhanced, and a long range correlation between dipole mo- 
ments is to be established as well, making the relaxation of photon-created clusters 
slower and slower. 

4. Summary 

In conclusion, based on a coupled photon-phonon model, we conduct a theoretical 
survey on the relaxation dynamics of nano-sized ferroelectric cluster in the paraelec- 
tric phase of BaTi03, which is experimentally observed as a real time correlation of 
speckle pattern between two soft X-ray laser pulses. The density matrix is calculated 
by a perturbative expansion up to the fourth order terms, so as to characterize the 
time dependence of scattering probability. We show that this time dependence is 
determined by the relaxation dynamics of photo-created ferroelectric cluster, which 
is manifested as a central peak in the phonon spectral function. Near the T c , cluster 
excitation is stable with a long relaxation time. While, at high temperature, it is 
suppressed by the thermal fluctuation, resulting in a short relaxation time. These 
features are in good agreement with Namikawa's experimental result. 
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